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We show that the coercive field in ferritin and ferrihydrite depends on the maximum magnetic 
field in a hysteresis loop and that coercivity and loop shifts depend both on the maximum and 
cooling fields. In the case of ferritin we show that the time dependence of the magnetization also 
depends on the maximum and previous cooling fields. This behavior is associated to changes in 
the intra-particle energy barriers imprinted by these fields. Accordingly, the dependence of the 
coercive and loop shift fields with the maximum field in ferritin and ferrihydrite can be described 
within the frame of a uniform-rotation model considering a dependence of the energy barrier with 
the maximum and the cooling fields. 

PACS numbers: 75.30.Cr, 75.50.Ee, 75.60.Ej, 75.50.Tt 



I. INTRODUCTION 



The magnetic properties of ferritin have been exten- 
sively studied in the last decade due to their puzzling 
features such as the existence of a maximum in the mag- 
netization derivative at zero field, [H,[3 a nonmonotonic 
field dependence of the magnetic viscosity, [lHH and a de- 
crease of the antiferromagnetic susceptibility with tem- 
perature below the Neel temperature when considered at 
low fields. [H [H Many of these studies were performed to 
enlighten the possible existence of quantum tunnelling 
in ferritin in the kelvin range. Less attention has been 
paid to the magnetic hysteresis, although it was termed 
"anomalous" in an early report due to the large coerciv- 
ity (~ 1800 Oe at 5 K), irreversibility up to relatively 
high fields (~ 35 kOe) and loop displacement found af- 
ter field cooling. Q Ferritin consists of a hollow spheri- 
cal shell composed of 24 protein subunits surrounding a 
ferrihydritc-like core. The diameter of the cavity is of 
the order of 7-8 nm and the average size of the core of 
horse spleen ferritin is 5 nm. Q The ferritin magnetic core 
orders antiferromagnetically and has both compensated 
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and uncompensated spins, resulting in a net magnetic 
moment of about a hundred of Bohr magnetons \ib per 
particle. 

Ferritin belongs to a group of nanoparticles where loop 
displacements are found but where strictly speaking there 
is neither a ferromagnetic (FM) material coupled to an 
antiferromagnetic (AF) one nor a cooling across a Neel 
transition temperature. The cooling is normally per- 
formed across the blocking temperature of the nanopar- 
ticles, i. e., across the temperature below which the mag- 
netic moment of the "average" nanoparticle is not able to 
fluctuate across the anisotropy barrier within the char- 
acteristic time of the measurement, which in ferritin and 
for dc measurements is about 20 K. In these systems, the 
origin of the loop shift is not clearly established. In AF 
nanoparticles the loop shift has been often interpreted 
as an exchange bias between the AF core and the un- 
compensated spins of the spin-glass surface. In the case 
of FM and fcrrimagnctic nanoparticles the loop shift is 
thought to arise due to exchange between the spin-glass 
spins and the FM core. [8, 9] Other studies attribute the 
loop shift to the existence of a minor hysteresis loop. [l(| 
In fact, as highlighted in Ref.@, Q, many of these sys- 
tems show high field irreversibility and non-saturating 
hysteresis loops, raising the question of the influence of 
the minor loops on the exchange field and, eventually, 
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the question of the existence of a real exchange bias. 

Here we report a study on the coercive field and loop 
shifts in ferritin at low temperature, obtained in mag- 
netization cycles recorded for different maximum fields 
up to 30 x 10 4 Oe, and after cooling under the influence 
of fields with different intensities. This study is also ex- 
tended to fcrrihydritc nanoparticlcs, which are similar 
to the ferritin magnetic cores, and complemented by the 
measurement of the time dependence of the magnetiza- 
tion near zero field. 



II. BACKGROUND 

Within the framework of the uniform rotation models 
developed by Stoner-Wohlfarth and Neel-Brown, [TTJ, [l2[ 
reversal of magnetization over an energy barrier E sepa- 
rating two minima is a coherent process, which can occur 
at T = K when the barrier is lowered down to zero by 
applying a magnetic field in the opposite direction of the 
particles magnetization. (l3| or it can occur by thermal 
activation when the thermal energy kgT becomes com- 
parable to E and thus the characteristic reversal time r 
becomes comparable to the characteristic measurement 
time r m ,[ll[. These quantities are related by the Arrhe- 
nius law 



r = tq exp 



E 



(1) 



where To is the inverse of an attempt frequency, supposed 
constant for simplicity in many situations. The energy 
barrier is field-dependent and it can be written as 



E{H) = E 1 - 



H 

#0 



(2) 



with a = 2 for systems with uniaxial anisotropy and easy 
axes parallel to the applied field, where E is the energy 
barrier at zero field and Hq is the switching field at zero 
temperature. In ferromagnetic materials Ho = 2K/Ms 
and Eq = KV, while in general 
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KV P , 



(3) 



where V is the particle volume, K is the anisotropy 
constant, Ms is the saturation magnetization and p 
an exponent equal to 1/2 in the case of antiferromag- 
netic ferrihydrite nanoparticles. [pi! In the framework 
of the Neel model and for a random distribution of 
anisotropy axes, a = 4/3 [Til ]. In general, simulations 
and experimental results show that a depends on the 
anisotropy, distribution of particle sizes and on intcrpar- 
ticle interactions. (lH Il7j 

According to Eq. magnetization reversal occurs 
at the coercive field He, when the energy barrier E{Hc) 
becomes small enough to be overcome at the given T and 
measuring time r 



H c 



1 



En 



1/a 



(4) 



where we have just re- written Eq. ([2]). The dependence 
of He on temperature, nanoparticle volume and charac- 
teristic measurement time can then be obtained by using 
Eq. (TD for E(Hc) in the previous equation, as shown in 
Ref . [lH Il9j . In particular, for the volume dependence of 
He at constant temperature one has 



H c (V)=H \l-(V B /V)W^ 



(5) 



where Vb is the blocking volume, i.e., the volume above 
which Eq cannot be crossed within r m in a system with 
anisotropy K and at a temperature T. A more refined 
expression for He (T) can be obtained b y co nsidering the 
temperature dependence of K and M,g.[20| 



III. EXPERIMENTAL 

Horse spleen ferritin samples used in these experiments 
were obtained from Sigma Chemical Company and pre- 
pared in powder samples by evaporation of the solvent 
at room temperature. 

The synthesis of the ferrihydrite nanoparticles in the 
organic-inorganic matrix (termed di-ureasil) has been de- 
scribed elsewhere. The particles are precipitated by 
thermal treatment at 80 °C, after the incorporation of 
iron nitrate in the matrix. The sample here studied has 
an iron concentration of 2.1 % in weight and the particles 
have a diameter distribution that can be described by a 
lognormal function 



f(D) 



1 



Ds D V2w 



exp 



{\og{D/n D )f 



2s D 



(0) 



with no = 4.7 ± 0.2 nm and standard deviation of the 
natural logarithm of the diameter = 0.43 ± 0.05. [22j 

For ferritin, magnetization was measured as a func- 
tion of field up to different maximum fields H max (in the 
0.5 x 10 4 to 30 x 10 4 Oe range) and after cooling from 100 
K down to low temperature (3.2 and 4.2 K) in zero-field 
cooling (ZFC) and field-cooling (FC) procedures using 
different cooling fields H coo i. These measurements were 
performed in a PPMS system (Quantum Design) with a 
vibrating sample magnetometer (VSM) option, and in a 
Bitter magnet with a VSM (HFML facility, Nijmcgen). 
In the latter, the modulus of the magnetization was mea- 
sured and magnetization curves were reconstructed by 
using the proper signal. Near zero this procedure is not 
perfect (since noise is always additive) and a small kink 
around zero field appears (see Fig. [1]). 

At the characteristic time of measurement, irreversibil- 
ity phenomena vanish for T > 40 K and magnetization 
curves taken at 3.2 K after cooling with H coo i = 2 x 10 4 
Oe from temperatures between 50 and 300 K are simi- 
lar. The magnetization curves are also independent on 
the cooling rate (cooling from 100 K) in the 0.5 to 5 
K/min range. In addition, magnetization was measured 
as a function of time during about 1000 s [M (t)] after 
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cooling from 100 K down to 4.5 K in ZFC and FC with 
H coo i = 0.5 x 10 4 Oe procedures. For each cooling pro- 
cedure and at low temperature (4.5 K), we have followed 
two different measurement protocols: i) applied different 
H m ax, removed the field down to a value near zero (50 
Oe) and measured M(t) and ii) applied different H max , 
then reversed the field to —H max , removed it down to 
—50 Oe and measured M(t). 

For the ferrihydritc nanoparticlcs grown in the organic- 
inorganic hybrid matrix, magnetization was measured as 
a function of field up to different maximum fields H max 
(in the 2 x 10 4 to 20 x 10 4 Oe range) and after cool- 
ing from 100 K down to 3.2 K in FC procedure using 
H C ooi — 2 x 10 4 Oe, in a Bitter magnet with an extrac- 
tion magnetometer (GHMFL facility, Grenoble). 



IV. RESULTS AND DISCUSSION 

A. Effect of H m ax and H coo i in the magnetization 
loops 

The magnetization loops of ferritin obtained at 
low temperature are analogous to those previously 
reported,]!], 0, HH with non-saturation, high-field irre- 
versibility and a maximum in the magnetization deriva- 
tive at zero field. The ZFC cycles are symmetric and in- 
creasingly broader as H max increases [Fig. HJa)]. When 
the sample is cooled in the presence of a field H coo i the 
magnetization curves are shifted in the H axis [Fig. QJb)] , 
being also increasingly broader as H max increases up to 
fields of the order of 10 x 10 4 Oe. Another interest- 
ing observation is that the differences in decrcasing-field 
branches of the loops obtained after a FC procedure for 
different H max are less significant than those obtained af- 
ter a ZFC procedure, while in increasing-field branches, 
the differences arc more significant after FC than after 
ZFC procedure. Similar results are found for ferrihydritc 
nanoparticles. 

The horizontal shift of the hysteresis loops in the 
cooling field direction is similar to that previously 
found in ferritin, ferrihvdrite [zij and other magnetic 
nanoparticles. [1] In the case of AF NiO nanoparticles, 
this loop shift was associated to surface anisotropy and 
multisublattice states, with the latter being associated to 
a variety of reversal paths. [2^] Surface anisotropy arises 
due to the breaking of the crystal-field symmetry at the 
boundary of the nanoparticle. Two models have been 
considered: one where the easy axis is transverse to the 
boundary and another where the local easy axis depends 
on the site "defect" (Neel surface anisotropy model). [2^ 
The loop shift in AF nanoparticlcs is also interpreted in 
terms of an exchange bias between antiferromagnetic and 
uncompensated moments, although no transition tem- 
perature is crossed but rather a blocking temperature. 

The effect of H max and H coo i on the field values at 
which magnetization crosses zero in the decreasing and 
increasing field branches of the hysteresis loop (termed 
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FIG. 1: Detail of magnetization loops for ferritin obtained at 
3.2 K (4.2 K in the case of Hmax = 30 x 10 4 Oe) after ZFC 
(a) and FC with H coo i = 2 x 10 4 Oe (b) measured for different 
H max - Hq and Hq correspond to the field values at which 
magnetization crosses zero in the decreasing and increasing 
field branches of the hysteresis loop, respectively. 



Hq and Hq , respectively) can be observed in Fig. [3] 
For Hcooi 7^ 0, Hq has a smaller variation with H max 
compared to that of Hq , while for H coo i = 0, Hq and Hq 
have symmetric variations. ^From Fig. [5J it is also clear 
that the effect of H max on Hq is more important than 
the effect of H coo i. In Fig. [3l we show the dependence of 
Hq, Hq, the coercive field He = {Hq — Hq)/2 and the 
loop shift Hs = — (Hq +Hq)/2 on the cooling field H coo i 
for the highest maximum applied field H max = 30 x 10 4 
Oe. With increasing H coo i, Hq increases approaching 
the ZFC value for H coo ; > 10 4 Oe, whereas \Hq\ values 
slightly decrease and arc always higher in modulus than 
the value for ZFC. Interestingly, the larger departures of 
Hq and H+ fr om the ZFC value occur for lower H coo i. 
As a result, Hq is almost independent on H coo \, being 
higher than the ZFC value, while the loop shift Hs has 
a small decrease with H rnn i. 
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FIG. 2: Fields at which magnetization crosses zero (Hq and 
Hg) for different cooling fields H coo i as a function of the maxi- 
mum field Hmax, for ferritin (a) and ferrihydrite nanoparticles 
(b). Lines represent fit to Eq. (I10I) . 



B. Effect of Hma X and H coo i in M(t) 



To have a better insight on the changes occurring near 
Hq and Hq induced by the application of different H max 
and Hcooi-, we have performed measurements of M(t) 
near zero field, after ZFC and FC under H coo i = 0.5 x 10 4 
Oe as described in Sec. IIIIl With these measurements, 
we aim at demonstrating that the changes in the hystere- 
sis loops are related to changes in the energy barriers to 
magnetization reversal induced by H max and H coo i- 

Before discussing the results of the magnetic relaxation 
measurements presented in Fig. EJc), we will start by an- 
alyzing the magnetic state of the samples attained after 
the protocols previous to the relaxation measurements. 
First, we plot in Fig. HJa) the values of magnetization 
obtained at H max (M@H max ) after ZFC and FC pro- 
cedures. As it can be observed, the magnetization val- 
ues obtained after a FC are always higher than those 
attained after a ZFC process. Second, the magnetization 
values obtained right after the field is decreased down 
to 50 Oe (-50 Oe) [M(t )} after FC and ZFC procedures 
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FIG. 3: Dependence of the fields at which magnetization 
crosses zero (Hg and Hq ) (a), He and Hs (b) on the cool- 
ing field Hcooi for a maximum field H max = 30 x 10 4 Oe, for 
ferritin. The ZFC values are shown in open symbols (which 
in the case of Hq and Hg corresponds to the same value). 
Lines connecting the FC values are eye-guides. 



are plotted in Fig. 0Jb). Comparison with results in 
the previous panel shows that the difference between the 
ZFC and FC values of M@H max (full squares) is always 
smaller than the differences between the ZFC and FC 
values of M(to) (open circles). In fact, at the highest 
applied H max = 5 x 10 4 Oe, FC and ZFC values of 
M@H max are identical, while the corresponding values 
for M(irj) become substantially different. The differences 
after FC and ZFC procedures that appear in Af (to) after 
decreasing the field down to ±50 Oe are an indication 
of the different energy barriers that each particle mag- 
netic moment has been able to cross. The higher M(to) 
values measured after FC hint at the appearance of en- 
ergy barriers imprinted after FC that are higher than af- 
ter the ZFC process. Moreover, the fact that M@H max 
values become the same after measuring at the highest 
H m ax = 5 x 10 4 Oe [see Fig. Ufa)] indicates that the dif- 
ferences in M(to) appearing after both procedures can- 
not be attributed to an increase of the net magnetic mo- 
ment of the individual particles induced at H max after 
the FC process. Finally, it is also interesting to note the 
constancy of M(to) observed for positive H max after FC 
[filled circles in Fig. Bfb)], which shows that the fraction 
of magnetic moments that reverse after reduction of the 
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field is almost independent of the maximum applied field 
and reinforces the two points commented previously. In 
contrast, when H max is applied in a direction opposite to 
Hcooi [points with negative abscissas in panels (a) and (b) 
of Fig. [4], there is a progressive increase of M(to) with 
increasing H max for both ZFC and FC procedures while 
M@H max values are essentially the same in the two cases. 
Again, this shows that the changes in M(to) cannot be 
attributed to an increase of the uncompensated moment 
of the particles but rather to the fact that, when arriving 
near zero field, different fractions of magnetic moments 
are able to cross the energy barriers at to depending on 
H max and H coo i, an indication that the effective energy 
barriers felt by the particles near zero field are modified 
by H max . In a simple picture, a negative H max has the 
effect of erasing the barriers imprinted by the positive 
H coo i, with the system being closer to the ZFC configu- 
ration as the intensity of H max increases. 

In what follows, we will analyze the results of the relax- 
ation measurements following the two above mentioned 
protocols. As previously found in ferritin[l|, M(t) dis- 
plays a quasi-linear dependence on ln(i) at intermediate 
times within the studied t range (up to 1000 s) and can 
be fitted to the following expression: 



M(t) ps Mq - Sln(t) 



(7) 



where Mq is related to the initial magnetization and 
S is the so-called coefficient of magnetic viscosity. Eq. 
([7| is particularly useful in situations where To is not 
known. (27| This equation can be derived from a general 
expression for the time dependence of the magnetization 
of an ensemble of nanoparticles with distribution of en- 
ergy barriers f{E), after field removal 



M(t) 



M (E)exp(-t/r(E))f(E)dE , (8) 



where Mq(E) is the initial magnetization of a particle 
with energy barrier E. It can also be shown that S is pro- 
portional to Ef(E) and, therefore, S is an appropriate 
quantity to observe changes in the energy barriers. [27i [28j 
However, the direction of this change is not directly given 
by S since, in principle, f(E) is a non monotonous func- 
tion. A distribution of energy barriers results directly 
from a distribution of volumes, according to the relation 
f(E) = g(V)(dV/dE), where g(V) is the volume distri- 
bution. Other sources of a distribution of energy barriers 
are a distribution of shapes, the existence of nanopar- 
ticles with the same size but with different degrees of 
crystallinity, different oxygen and water content. 

This equation assumes that the magnetization decay 
of a nanoparticle ensemble is due to the switching of the 
nanoparticles magnetic moments as a consequence of en- 
ergy barrier crossing when, for a given T and H, the 
Arrhenius relaxation time r is of the order of the mea- 
surement time r m . 

The H max dependence of the viscosity coefficient S as 
obtained from fits of the linear part of the relaxation 



curves to Eq. [7] is reported in Fig. |4jc). We observe 
that for positive H max and for the FC case, S remains 
essentially constant with increasing H max . This means 
that H C ooi imprints energy barriers for reversal into the 
FC direction that are not substantially changed by a pos- 
itive applied H max . However, in the ZFC case, the re- 
laxation rate increases with H max , showing that H max 
changes the energy barriers in this case. For negative 
H m ax (applied contrary to the H coo i direction), however, 
the energy barriers are shifted. In summary, the general 
behavior of S is similar to that of M(to), reinforcing the 
interpretation of the effects of H coo i and H max in terms 
of energy barriers. Since in ferritin (29. 1301 and in the 
fcrrihydritc nanoparticles here studied 14| intcrparticle 
interactions are negligible, the proposed changes in the 
energy barriers are most probably associated to intra- 
particlc phenomena, as discussed below. 



C. Effect of H max on the energy barriers 

In Fig. [SJ we present the dependence of He and 
Hs with H max for hysteresis loops measured after ZFC 
(squares) and after cooling in different different H coo i- 
The first point to notice is that, in the FC case, the loop 
shift Hs first rapidly decreases for low H max but, for 
fields higher than 10 x 10 4 Oe, it saturates to a value 
around 500 Oe even for extraordinarily high values of 



Hn 



30 x 10 Oe. This behavior is somewhat unex- 



pected since, for the usually observed loop shifts due to 
minor loops, the shift tends to zero for sufficiently high 
fields. 

Secondly, the constancy of He obtained at high fields 
in the FC case is an indication of the existence of high 
energy barriers imprinted by H coo i that cannot be sur- 
mounted even by applying a H max of 30 x 10 4 Oe field in 
the direction opposite to H coo i- Both observations (to- 
gether with the Hq variation already presented in Fig. 
[3]) indicate an initial evolution of the minor loops due 
to crossing of smallest energy barriers. After this first 
stage, the variations are smoothed by the higher energy 
barriers imprinted by H coo i- According to the behavior 
of Hq, Hq and M(to), it is clear that H coo i increases the 
energy barriers in the field direction (Eq) and decreases 
barriers in the opposite direction (Eq). Since the high- 
est barriers are not overcome, the symmetric situation 
Eq = Eq cannot be recovered and thus Hs is always 
different from zero. 

A quantitative description of these phenomena can be 
given within the framework of the uniform rotation model 
in terms of the influence of H max on the energy barriers 
near zero field Eq. First, let us notice that the depen- 
dence of He on H max is similar to the dependence of 
He on V usually found in nanoparticle systems (see Eq. 

Taking into account this resemblance, and with 
the aim to propose an expression for He(H max ) which 
properly describes the measured data, we will assume 
that H max influences the zero field energy barriers in a 
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FIG. 4: (Colour online) (a) Magnetization measured at 
Hmax{—H ma x) before removing the field down to 50(— 50) Oe 
(M@Hmax) obtained after ZFC and FC with H cool = 0.5 x 10 4 
Oe, and difference between M@H ma x after FC and ZFC 
AM@H ma x', (b) first value of the remanent magnetization 
measured at 50(-50) Oe [M(t )] obtained after ZFC and FC 
with Hcooi = 0.5 x 10 4 Oe and difference between the FC and 
ZFC values; (c) magnetic viscosity S measured at 50(— 50) 
Oe after a previous H ma x{—H ma x)- Measurements were per- 
formed in ferritin. 



cles. Therefore, the influence of H max can be thought 
mostly as an effect on K associated to an increase of the 
local (intra-particles) energy barriers. Since these are 
macroscopically average measurements, it is difficult to 
access the "microscopic" origin for this effect on K . Any- 
way, this can be understood considering that the system 
has multiple configurations with associated energies such 
that H m ax and H coo i selects or imprints a set of these 
configurations restricting the relaxation of the moments. 

Based on Eq. Q the relation between Hq and H max 
can be expressed as 



H+{H max ) = H+{oo) 



1 - 



Hi 
H m 



0+' 



(10) 



where (3 + = a/7, Hq(oo) is redefined as Hq for infinite 
Hmax and Hg is defined as the field at which Hq is zero. 
The relation between Eq and Hg is that expressed in 
Eq. © since H^ is a particular case of an H max . When 
comparing to experimental data, care must be taken and 
the FC and ZFC cases must be distinguished. In the FC 
case, Hq can be negative and Hg is the field at which 
Hq crosses zero (Fig. [5]). In the ZFC case, Hg is defined 
as the field at which the behavior of Hq at high H max 
extrapolates to zero, since in practice the experimental 
Hq values are not zero at H max ^ H^. 

Eq. (|TU|) can successfully describe the Hq (H max ) 
data shown in Fig. [2j with H£(oo) = 0.23 x 10 4 Oe, 
H+ = -1.1 x 10 4 Oe and /3+ = 0.8 for the FC data. 
In the case of ferrihydrite Hq(oo) — 0.20 x 10 4 Oe, 
H+ = -2.0 x 10 4 Oe, p + = 1.5. Considering a = 4/3, 
7 rj 1 in ferritin and 7 ~ 2 in ferrihydrite, i. e. approx- 
imately a linear and quadratic relation between Eq and 
Hmax- The differences here found for 7 are associated 
to the fact that in ferrihydrite a smaller H max is enough 
to approach Hq to saturation. Again, the "microscopic" 
origin for this mechanism is not clear. 

For ZFC data, both HQ{H max ) and HQ~(H m ax) are 
well described by Eq. (JTDJ) for |i? maa: | > \Hg\, while for 
FC Hq (H max ) is approximately constant. Accordingly, 
the generalization to the H max dependence of the coer- 
cive field for ZFC and FC procedure (Hczfc and Hcfc, 
respectively) and of Hs is straightforward 



way similar to the particle volume (Eq. [3]) 

£0 « H max (9) 

where 7 is a power law exponent that controls the way 
He approaches its limiting value for high H ma x, such 
that a higher 7 is associated to a faster approach to satu- 
ration. This power law dependence condensates different 
possible mechanisms for the influence of H max on intrin- 
sic energy barriers: either a change in the exponent p 
in Eq. Q or an irreversible increase of the anisotropy 
constant K. The first possibility seems to be ruled out 
since, as reported in the previous section, H max does not 
seem to affect the net magnetic moment of the parti- 



Hc ZFC (Hmax) = Hq (oo) 



H 



CFC 



l -H 

2 



II 




and 



H s (H max ) = --Hq --H+(oo) 



1 - 



H 



H„ 



(12) 

where H is constant. 

As expected from the agreement between the 
H^Hmax) data and Eq. (JTUJ), Eq. (fTTj) can also be 
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successfully used to describe the HczFc{H max ) and 
H C Fc(H max ) data in the 0.5 x 10 4 < H max < 30 x 10 4 
Oe range, as shown in Fig. EJa), with Hq = 0.24 x 10 4 Oe 
and H+ = 0.9 x 10 4 Oe for the ZFC data, and with the 
previous Hq (H max ) parameters and Hq = —0.297 x 10 4 
Oe for the FC data. For low fields {H max < 0.5 x 10 4 
Oe) the fit deviates from the experimental data, while 
the fit considering only the high field data extrapolates to 
Hq = 0. At Hg, the experimental He is still about 18% 
of its saturation value, approaching zero for H max = 0. 
The differences between the Hcz fc and He fc data are 
probably of the order of data error and, thus, the dif- 
ferences between ZFC and FC fitted parameters are also 
within the error bars. 

The value of He at H max = 5 x 10 4 Oe and 3.2 and 
4.2 K is of the order of that previously found for ferritin 
at 5 K (~ 1700 Oe).[|| The slightly higher value that was 
found (2200 Oe) is probably due to the lower tempera- 
ture. Other factors affecting He that may contribute to 
this difference are the field sweeping rate and character- 
istic time of measurement. 

As in the case of He , the dependence of H$ with H max 
can be in fact described by Eq. (|T2"T) , with the parameters 
obtained for Hq and a constant Hq = —0.297 x 10 4 Oe, 
which gives an extrapolated H$ (oo) = 335 Oe. 




O.OO r 

5 10 15 20 25 30 



FIG. 5: Dependence of the coercive field He (a) and loop shift 
Hs (b) on the maximum field H ma x in ferritin, for different 
cooling fields H coo i. Lines in panels (a) and (b) represent fits 
to Eq. (|lip and Eq. (|f 21) . respectively. Inset shows zoom over 
the low H m ax region. 

V. CONCLUSIONS 

Coercivity and loop shifts in nanoparticles are dynam- 
ical phenomena, which depend on temperature, charac- 



teristic time of measurement and number of cycles, for 
instance. In ferritin, we have shown that coercivity and 
loop shifts depend also on the cooling field and on the 
maximum field used, for fields higher than those normally 
used. The dependence of coercivity and loop shifts with 
the maximum field can be described in terms of changes 
in the anisotropy energy barrier near zero field induced 
by the maximum field, and quantitatively described by 
a modified Neel-Brown model here proposed. Qualita- 
tively, field cooling imprints energy barriers, such that 
the energy barriers near zero in the descending and as- 
cending branches of the magnetization cycle are higher 
and lower than in the ZFC case, respectively. This dif- 
ference is attenuated (but not erased) by increasing the 
maximum field in the opposite direction of the cooling 
field. Accordingly, the loop shift decreases with the max- 
imum field but it is not zero up to the highest field used 
(30 x 10 4 Oe), showing that the barriers imprinted by 
field cooling cannot be overcome by these high fields. 

The experimental observations and subsequent analy- 
sis presented in this article have evidenced the imprinting 
of high energy barriers through an effective anisotropy 
induced by the applied protocols. This gives rise to 
the high irreversibility and minor loop effects similar to 
those observed in spin-glasses and diluted antiferromag- 
nets, where this phenomenology is ascribed to dilution 
and the antiferromagnetic character of the interactions 
and not to frustration. [3ll |32| 
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